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SUMMARY 


We used genetic data from nine populations of Grindelia howellii, one population of 


Grindelia squarrosa, and one putative hybrid population to examine genetic diversity and 
variability. The general results are provided in the table below: 























Analysis Question Results 
: Are the markers Data quality is good - no evidence of missing 
Data quality : : : : 
appropriate? genetic data or correlations between loci 
: : pte BOP ulauicns Measures of inbreeding are not particularly high; 
Inbreeding inbred and in poor j : : 
: SO-148 has the highest measure of inbreeding. 
genetic health? 
, Is there high genetic Genetic diversity is not low in the G. howellii 
Genetic . eae : ; 
ae dk diversity in the populations, but we don't have enough G. squarrosa 
diversity : : eer ; 
populations? populations to make a determination in that species. 
The combined populations, including both species, 
oer re are distinct from each other. G. howellii populations 
Genetic 6 San from the Seely Lake area (Montana populations) are 
cate genetically similar to nee : 
variability very similar to each other. Most of the genetic 
each other? pati, Mt ye tee ew 
variability is within individuals, not among 
populations or individuals 
Genetic Are geographically There is a correlation between genetic and 
versus distant populations geographic distance across all study populations. 
geographic genetically distant as | Within G. howellii Montana populations, there is no 
distance well? correlation. 
. ee oe nan We don't have enough samples of G. squarrosa and 
Hynes na walls utative hybrids to answer this question 
hybridizing? P y 4 ' 





We used the data to answer the following questions. 


(1) Are the G. howellii populations’ sizes, genetic diversity, and variability ideal to preserve 


genetic integrity? 


The number of total individuals in the sampled populations varies greatly (30 — 8,800), 
but genetic diversity is not strikingly different across populations. If we accept the general 
northern, central, and eastern groupings of the Seely Lake (Montana) populations as ‘meta- 
populations’, the northern meta-population has lower genetic diversity than the central and 
eastern meta-populations. If the northern populations are large, the relatively low genetic 
diversity could be of concern. 


(2) How do the genetics of G. howellii interact with conservation efforts, restoration, and 
management for this species? 


Our results do not indicate any immediate threat to the sampled populations in terms of 
low genetic diversity, and indeed the Montana populations seem to be acting as a large meta- 
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population, with few genetic differences between populations. In particular, SO-143 is found 
along a road that may be developed in the future. Looking at its genetic diversity and the 
STRUCTURE results, it is not distinctly different from SO-148 and SO-82/83. However, it does 
have a private allele found only in that population. If the population is destroyed, it may be worth 
transplanting individuals and/or collecting seeds to add to the SO-148 and SO-82/83 populations. 


(3) How closely related are G. howellii and G. squarrosa? 

Species identification within Grindelia is hard due to the wide range of morphological 
variation within and among species, and previous genetic studies have found that ‘species’ aren’t 
always unique genetic units. The putative hybrid population in this study, SO-171, is quite 
distinct from the other G. howellii Montana populations. However, if it is a hybrid, we’d expect 
it to share genetics with both the Blue Mountain population of G. squarrosa and Montana 
populations of G. howellii. However, Blue Mountain, though different from the Montana 
populations, isn’t strikingly so — as a different species, we’d expect it to be quite distinct. The 
SO-171 population could be a hybrid, or it could be that geographic distance from the other 
Montana populations has led to genetic differentiation. The Clear Creek population is also quite 
distinct. This could be because of geographic distance or because it is not G. howellii. At this 
point, we recommend collecting from more populations of G. squarrosa and putative hybrid 
populations. 


(4) What is the relationship between geographic range size and levels of genetic diversity in 
lineages? 

Without more G. squarrosa populations, we can’t relate population genetic levels to 
range sizes. The single included putative G. squarrosa population doesn’t have markedly 
different levels of genetic diversity compared to the G. howellii populations. 
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INTRODUCTION 

Grindelia howellii (Asteraceae; Howell’s Gumweed) is a restricted species endemic to 
Montana and Idaho (USDA National Resources Conservation Service, 2017). There is debate 
over classification of the species; G. howellii has been grouped with Grindelia paysonorum 
(Bartoli and Tortosa, 2012) and considered similar to Grindelia nana (Strother and Wetter, 
2006). However, as Strother and Wetter (2006) note, G. howellii has distinctive glands on the 
stem and leaves, making it sticky to the touch. The most recent molecular phylogeny of the 
genus confirms its place in the “Pacific” Clade with other species in the region (Moore et al., 
2012). However, this genetic study only tested one herbarium specimen (accession) for G. 
howellii, and other species with multiple accessions per species were polyphyletic. These results 
point to a complicated relationship between taxonomy and evolutionary differences in the genus. 

Beyond these taxonomic issues, G. howellii is a considered a sensitive species in 
Montana by the US Forest Service. It is found in open habitats such as roadsides and open 
wetlands (Lesica, Lavin, and Stickney, 2012), and could be negatively impacted by road 
management or development. Also, given its restricted range, populations may have low genetic 
diversity and/or require protection in order to perpetuate the species. Despite local abundance 
(Lesica, Yurkewycz, and Crone, 2006), restricted endemics can be limited by dispersal distances 
and low genetic diversity. Their limited range and habitat preference can make them more 
susceptible to fragmentation and loss of genetic diversity (Wagenius, Lonsdorf, and Neuhauser, 
2007). However, endemism and restricted range does not always indicate low genetic diversity 
or high inbreeding (Williams et al., 2016). In addition to occupying a small range, G. howellii 
may hybridize with other species in the genus, such as Grindelia squarrosa (Lesica, Lavin, and 
Stickney, 2012), which is widely distributed across the United States. Grindelia squarrosa has 
reportedly hybridized with at least one other species, G. nana (Moore et al., 2012). Hybrid 
populations may exist where species overlap in their ranges, and produce morphologically 
intermediate individuals. 

Here, we used microsatellite loci and samples from G. howellii and G. squarrosa 
populations to examine genetic diversity and possible hybridization. The main questions to be 
answered in this study were: 1) Are the G. howellii populations’ sizes, genetic diversity, and 
variability ideal to preserve genetic integrity? 2) How do the genetics of G. howellii interact with 
conservation efforts, restoration, and management for this species? 3) How closely related are G. 
howellii and G. squarrosa? 4) What is the relationship between geographic range size and levels 
of genetic diversity in G. howellii and G. squarrosa? 
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MATERIALS AND METHODS 


Sampling 

Leaf samples were collected from nine populations of G. howellii, one population of G. 
squarrosa, and one population of a putative G. squarrosa X G. howellii hybrid' (Fig. 1; Table 1). 
The Clear Creek population of G. howellii occurs in Idaho while the remaining populations occur 
in Montana. At each population, 3-5 leaves per individual were collected from 16 to 33 
individual plants; thus, a total of 280 plants were sampled. The only exception to this occurred 
at the G. squarrosa “Blue Mountain” population, where 3-5 leaves were collected from 7 
individual plants. We ran analyses on three different datasets: all samples from all 11 
populations, just G. howellii (excluding G. squarrosa and the putative hybrid population), and 
just the Montana populations of G. howellii (which excluded the Clear Creek, Idaho population 
of G. howellii). 


Microsatellite Loci 

DNA was extracted from fresh leaf tissue using a DNEasy 96 Plant Kit (Qiagen Valencia, 
California, USA) and quantified using a ThermoFisher Quant-iT PicoGreen dsDNA Assay Kit 
visualized on a Molecular Devices Gemini XPS plate fluorometer. The DNA samples were 
genotyped using 6 microsatellite loci (Table 2). Loci were originally developed for Grindelia 
species (Abercrombie et al., 2009), and have been used successfully in other studies of Grindelia 
(Moore, Moore, and Baldwin, 2014). 

The loci were amplified with 10 ul PCR reactions (1.5U HotStar Taq, 0.2 mg/ml BSA, 
0.4 uM forward primer, 0.4 uM reverse primer, | ul of 10x buffer with 15 mM MgCl, 125 uM 
dNTPs, 10 ng DNA; 95°C for 5 min, [94°C for 30 sec, 55°C for 45 sec (-0.5°C per cycle), 72°C 
for 45 sec, repeat 20x] [94°C for 30 sec, 45°C for 45 sec, 72°C for 45 sec, repeat 10x], 72°C for 
10 min). The samples were run on an ABI3130x1 with ROX as an internal standard, with the 
GRINO26 amplification product diluted 1:30. The loci were scored by hand using GeneMarker 
v1.95 (SoftGenetics LLC, 2001-2010). 


Analyses 

We checked for null alleles using Micro-Checker with 1000 replicates (Van Oosterhout et 
al., 2004) and examined linkage disequilibrium (LD) among loci using Genepop (Raymond and 
Rousset, 1995; Rousset, 2008) with the following Markov chain parameters: 100 batches, 1000 
iterations per batch, and dememorisation of 1000. We also used Genepop to determine Fis values 
and their significance with default options. 

We used GenAIEx v6.503 (Peakall and Smouse, 2006; Peakall and Smouse, 2012) to 
calculate diversity statistics (number of alleles, effective number of alleles, observed and 
expected heterozygosity, private alleles, Fis, Fst, AMOVA), visualize Nei’s genetic distance 


1 The Montana Natural Heritage Program botanist proposed the idea that this population could be a hybrid 
because morphologically plants exhibit leaves that barely clasp the stem and have a lower density of stipitate- 
glandular hairs on the stem making them less sticky. 
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(Nei, 1972) using Principal Coordinates Analysis (PCoA) ordination, and run Mantel tests of 
geographic and Nei’s genetic distance with 999 permutations. 

We used the program STRUCTURE (Pritchard, Stephens, and Donnelly, 2000) to look 
for patterns of genetic structure. We used an admixture model with 1,000,000 generations with 
100,000 generations burn-in, with 5 iterations of K from one to the number of sampled 
populations — 1, which depended on the dataset. We used the method described in Evanno et al. 
(2005) as implemented in STRUCTURE HARVESTER (Earl and vonHoldt, 2012) to determine 
the optimal value of K, or number of genetic clusters in the data. 


RESULTS 


Locus Amplification 

Two loci used, GRIN24 and GRIN113, had evidence for null loci. For a large number of 
individuals (99 of 280 plants for GRIN24, and 20 of 280 plants for GRIN113), these loci did not 
amplify across multiple replicate PCR and ABI runs. For null loci, the entire locus fails to 
amplify, which leads to missing data. Null loci can result from a mutation in the primer 
sequence, which can be confirmed by sequencing the locus and primer region. To deal with the 
missing data in this study, we inferred loci frequencies when possible in GenAIEx. 

Locus GRIN68 was monomorphic for all individual plants except for one individual in 
population SO-171 from the putative hybrid population. This means that analyses that include 
only G. howellii populations effectively used five loci instead of six. 

Data Quality 

We also tested for null alleles, which are different from null loci. In the case of null 
alleles, a particular allele fails to amplify, which can lead to potentially high estimates of 
homozygosity. MicroChecker found no evidence for null alleles in this dataset. Genepop also 
found no evidence for linkage disequilibrium (p > 0.08); however, the relationship between 
GRIN68 and GRIN113, GRIN24, and GRIN45 could not be evaluated. This is likely because of 
the lack of variation in GRIN68. 


Population Genetic Diversity 

A combination of different measures can be used to evaluate population genetic diversity. 
Populations with high genetic diversity generally have low inbreeding coefficients, a high 
number of effective alleles, and similar measures of observed and expected heterozygosity. For 
all populations of Grindelia in this study, the inbreeding coefficient Fis (low values indicate less 
inbreeding) is quite variable (Table 3), with possible inbreeding at SO-148. The overall estimate 
of Fis across populations is 0.211. The number of alleles (Na) is similar across populations, 
which in turn is similar to the effective number of alleles (Ne). There is a similar pattern with 
observed heterozygosity (Ho) and expected heterozygosity (He), with slightly lower values of 
observed versus expected heterozygosity. Although individual loci within populations were 
significantly different from Hardy-Weinberg equilibrium, no population had all loci significantly 
different, and no locus was consistently significantly different across all populations. However, 
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GRIN113 was often significantly different. Sampling additional populations would allow us to 
evaluate if values of genetic diversity in G. howellii populations are significantly different from 
G. squarrosa. However, this analysis is not possible given that only one G. squarrosa population 
was sampled. 

One measure of population differentiation is Fst, which measures population 
differentiation and ranges from zero to one, with higher values indicating high population 
differentiation (Table 4). For all populations, Fst = 0.37 which indicates some differentiation. 
For the G. howellii populations, Fst = 0.28, which may reflect the genetic difference between 
Clear Creek and the SO populations. For the Montana populations, Fst = 0.13, which indicates 
little differentiation across the SO populations. These Fst measures mirror the results of the 
AMOVA tests. For all datasets, most variation was found within individuals. However, the 
Montana G. howellii populations had the least amount of variation among populations (13%) 
compared to all G. howellii populations (28%) and all populations (37%). In addition, private 
alleles (alleles found in only one population) can be an indication of population differentiation 
(Table 2). The Clear Creek population had the most private alleles, with 7, followed by the Blue 
Mountain population with 4. 


Population Structure 

The PCoA ordination (Fig. 2a) shows that the putative hybrid population SO-171, Clear 
Creek, and Blue Mountain populations are genetically different from the Montana G. howellii 
populations, and seems to correspond with geographic distance. Looking across the Montana G. 
howellii populations (Fig. 2b) there is loose cluster of the northern populations (SO-1, SO-172, 
and SO-52) and eastern populations (SO-11, SO-82/83, SO-143, and SO-148) with SO-115 
located in a somewhat different position. This is supported by the Mantel test results. Across all 
populations, there was a significant correlation between geographic and genetic distance (p = 
0.04, R’ = 0.9). A significant correlation between geographic and genetic distance is also found 
across all G. howellii populations (p = 0.04, R” = 0.43). However, no significant correlation 
between geographic and genetic distance is found for the Montana G. howellii populations alone 
(p = 0.4, R? = 0.0004). 

The STRUCTURE and STRUCTURE Harvester results show clear differences between 
SO-171, Clear Creek, and the rest of the populations, with K = 4 as the optimal number of 
clusters (Fig. 3). For all G. howellii, the most likely number of clusters is K = 2, with clear 
separation between Clear Creek and the SO populations. For just the Montana G. howellii 
populations, the most likely number of clusters is K = 3, which shows the distinction between 
northern populations (SO-52, SO-4, and SO-1), eastern populations (SO-82/83, SO-143, and SO- 
148), and the two central populations (SO-11 and SO-115). 


DISCUSSION 
In general, we found little evidence of low genetic diversity in populations of Grindelia 
howellii. Without more samples of Grindelia squarrosa, however, we are unable to compare 
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genetic diversity of G. howellii to a more wide-spread species or definitively identify a putative 
hybrid population. We discuss the specific project questions below. 


(1) G. howellii populations’ genetic integrity and variability 

Our data shows that the Clear Creek population of G. howellii in Idaho is distinct from 
the SO populations from Montana. It has a large number of private alleles, and consistently 
separates out from the other populations in the PCoA ordination and the STRUCTURE analysis. 
We suggest collecting more populations from Idaho and/or verifying the identification of the 
Clear Creek, Idaho population using Flora of North America and Bartoli and Tortosa (2012) 
botanical keys. 

Concentrating on the Montana G. howellii populations, our measures of diversity 
including number of alleles, number of effective alleles, and observed and expected 
heterozygosity are not particularly low. In terms of variability and population structure, we found 
that geographic distance between the Montana populations was not correlated with genetic 
distance, and Fst is low, indicating that these populations are genetically similar and likely 
exchanging genes though pollination or seed dispersal. We do see indications of some 
differences between populations using STUCTURE, with three general ‘meta-populations’. The 
northern meta-population has relatively lower genetic diversity that the central and eastern meta- 
populations. If the northern populations are large in terms of absolute number of individuals, the 
relatively low genetic diversity could be of concern. 

Inbreeding doesn’t appear to be a concern in most of the sampled populations. The 
exception is population SO-148, which has the highest measure of inbreeding and a low 
observed heterozygosity. Inbreeding may be low overall in G. howellii due to self- 
incompatibility which is widespread in the family Asteraceae (Charlesworth, 1985). Although a 
barrier to inbreeding, if self-incompatibility is present in the genus or species it could limit the 
number of potential mates in small populations. 


(2) G. howellii genetics as it relates to conservation efforts, restoration, and management 

The historic habitat of G. howellii is likely shallow wetlands and moist meadows (Lesica, 
Lavin, and Stickney, 2012), but it also currently favors disturbed habitats such as roadsides. Our 
results do not indicate any immediate threat to the sampled populations in terms of low genetic 
diversity, and indeed the Montana G. howellii populations seem to be acting as a large meta- 
population, with few genetic differences between populations. The sampled populations are 
found within a larger matrix of populations in the Seely Lake and Ovando valleys. 

In particular, SO-143 is found along a road that is proposed for re-construction. Looking 
at its genetic diversity and the STRUCTURE results, it is not distinctly different from SO-148 
and SO-82/83. However, it does have a private allele found in no other populations. It is a large 
population (estimated 8,800 individuals) which may harbor additional genetic diversity. If the 
population is to be destroyed, it may be worth transplanting individuals and/or collecting seeds to 
add to the SO-148 and SO-82/83 populations. 
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(3) Relationship between G. howellii and G. squarrosa 

Species identification within Grindelia is difficult due to the wide range of morphological 
variation within and among species, and previous genetic studies have found that ‘species’ aren’t 
always unique genetic units (Moore et al., 2012; Moore, Moore, and Baldwin, 2014). Moore et 
al. (2012) found that individuals from the same Grindelia ‘species’ did not group together using 
a phylogenetic approach. Previous work by Moore, Moore, and Baldwin (2014) using the same 
microsatellites in California populations found that populations segregated by ecotype, not 
ploidy level. In general, however, Moore et al. (2012) found that G. howellii and G. squarrosa 
were in different clades (genetic groups), with G. howellii in the Pacific Clade and G. squarrosa 
in a widespread Eastern United States Clade. 

However, there could be hybridization between the two species. The main barrier to 
determining the relationship between G. howellii and G. squarrosa and the existence of hybrid 
populations is that only one population of G. squarrosa and one population of a putative hybrid 
were included in this genetic study. The putative hybrid population in this study, SO-171, is quite 
distinct from the other Montana populations. However, if it is a hybrid, we’d expect it to share 
genetics with both the Blue Mountain population of G. squarrosa and the Montana G. howellii 
populations. However, Blue Mountain, though different from the Montana populations, isn’t 
strikingly so, and as a different species, we’d expect it to be quite distinct. The SO-171 
population could be a hybrid, or it could be that geographic distance from the other Montana 
populations has led to genetic differentiation. Also, the identification of G. squarrosa from this 
population is in doubt as both G. squarrosa and G. nana occur in the county (A. Pipp, pers. 
comm.). Future field work will confirm the identification of this population. 

Similarly, the Clear Creek population is also quite distinct. This could be because of 
geographic distance or because it is not G. howellii. We recommend collecting from more 
populations of G. squarrosa and putative hybrid populations. However, the SO-populations 
excluding SO-171 appear to be a cohesive genetic group with morphology consistent with G. 
howellii. 


(4) Geographic range size and levels of genetic diversity 

Without more G. sguarrosa populations sampled, we can’t relate population genetic 
levels to range sizes. The single included putative G. squarrosa population doesn’t have 
markedly different levels of genetic diversity compared to G. howellii populations. Future 
collections should be made of specimens with clear G. squarrosa characters and of specimens 
with clear G. nana characters, especially given that these species may hybridize in parts of 
Montana. 


CONCLUSIONS 

Grindelia howellii populations are not genetically depauperate compared to a single 
sampled population of a wide-spread congener species, G. squarrosa. The most genetically 
distinct G. howellii population in the study is from Clear Creek. The putative hybrid population 
SO-171 and the Blue Mountain (G. squarrosa) are also genetically distinct. The Montana G. 
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howellii can be divided into three meta-populations, north, central, and east, but appear to 
exchange alleles. There is currently no evidence for hybridization, but additional populations of 
G. squarrosa, G. nana, and putative hybrids are needed to fully test that hypothesis. 
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Table 1 — Eleven Grindelia populations were sampled for genetic analysis. 













































































‘ Estimated | Estimated # 
: Population : : : # of 
Species 1 Latitude, Longitude population of 
Name is date samples 
area (acres) | individuals 
G. howellii SO-1 47.438060, -113.645650 4 >300 23 
G. howellii SO-11 47.153300, -113.461790 13 >1340 31 
G. howellii SO-52 47.469470, -113.632820 <1 30 16 
G. howellii SO-82/83 47.073917, -113.219754 3 >500 29 
G. howellii SO-115 47.092840, -113.420140 <0.5 89 33 
G. howellii SO-143 47.110151, -113.184390 61 8800 29 
G. howellii SO-148 47.115588, -113.218901 47 1530 29 
G. howellii SO-1722 47.437780, -113.657373 <2 >100 26 
G. howellii x SO-171 46.502610, -112.531100 3 +300 7 
squarrosa 
G. howellii Clear Creek, 46.080200, -115.694900 NA NA 30 
Idaho 
G. Squarrosa Tee 46.825800, -114.099500 NA NA 7 
1 - All populations are in Montana unless noted 
2 - Originally labelled incorrectly in the field as coming from SO-4. It was collected from SO-172. 
3 - Averaged location for three sub-locations: (46.8288, -114.1), (46.8237, -114.102), (46.8249, -114.09) 
Table 2 — Microsatellite primers used in this study. 
Locus F Primer R Primer Boman Daner: || eoize 
dye Range 
TCGACTCGAATTCTC | TACATCGCCCAAATC 218- 
Ce ATGTGTC AATCA Nee 238 
CCATAGGCGGTATAC | CCCAATCATCCAACG 211- 
Pena AATAGCC AAATC aM 246 
GCTTCAAACTCAAGC | CACAAAAGCATGTGC 214- 
CRINGE AAGCA ATCAA ee 226 
TTGACTATCATTTGA | AGAGAAACATCCGTG 399- 
GRINO45 | ACCATCCA CGTGT NEP 419 
TGGGCTGACTCCTGC | TGCAGATTCCGAAAA 352- 
GRINS | TTAAT CACAA is 354 
CGTGGTCTAAACGGG | CCATTCCTCATCGAAC 454- 
Sate GTATC TCTTG anes 481 




















Note: All loci from Abercrombie et al. 2009, see text for annealing temperatures 
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Table 3 — Diversity statistics for eleven populations of Grindelia. 













































































Private Private Pavel 
Species Eepuanon Fis | Na | Ne | Ho | He Hardy Wandelt s alleles - all alleles - G.|__ Onin Notes 
Name Equilibrium ; .. |populations of 
populations | howellii in 
G. howellii 
G. howeltti | SO-1 | 0.297 |1.83|1.49]0.15 0.21] , Locus 11 1 1 1 3 of 6 loci are 
significantly different monomorphic 
G. howell’ | SO-11 | 0.219 |3.17]1.75|0.23/0.30]_. LOC 24 and 113 2 2 2 POL aS 
significantly different monomorphic 
G. howeltii | $O-52 |-0.146|1.50]1.26|0.17/0.14| Not significant 0 0 0 Pere tog! ae 
monomorphic 
G. howellit | SO-82/83 | 0.163 |2.67] 1.7 |0.26|0.31|_, LOC 24 and 113 1 1 1 ea 
significantly different monomorphic 
G. howeltit | SO-115 | 0.159 |2.00]1.43}0.20]0.24| Not significantly 0 0 0 ore eae 
different monomorphic 
G. howeltii | $0-143. | 0.240 |2.33]1.61|0.19|0.24| , _. Locus 113 1 1 1 2 of 6 loci are 
significantly different monomorphic 
G. howellii | $0-148 | 0.374 |2.83|1.89|0.22|0.34| , Locus 113 0 0 1 1 of 6 loci are 
significantly different monomorphic 
G. howelt | $O-172. | 0.161 | 1.67] 1.55]0.21 0.25] ___ Locus 113 0 0 0 Danan 
significantly different monomorphic 
© PONEEN A) 5694908) (1.83 | 1242/0160 24|'o See 1 NA NA one gee 
squarrosa significantly different monomorphic 
Clear : ; 
G. howeltii | Creek | 0.205 |2.83|1.58|0.26]0.32| , 10°! 26 and 45 7 7 a coke le ae 
daha significantly different monomorphic 
C spiaitrosal O° = 6495 |2'93 149 1007|0:90|' Novsemueantly 4 NA NA 1 of 6 loci are 
Mountain different monomorphic 
Note: Na = number of alleles, Ne = number of effective alleles, Ho = observed heterozygosity, He = expected heterozygosity 
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Table 4 — Measures of Fst and AMOVA results. 
































AMOV _| Within Among Among 
Dataset Fsr_| A results | individuals | populations individuals 
All populations ae p<0.01 50% 37% 13% 
G. howellii ' p<0.01 57% 28% 15% 
Montana G. howellii "5 p<0.01 69% 13% 18% 
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(b) 


Figure | — (a) Known Grindelia howellii populations in the Seely Lake Area (SO populations). 
Sampled populations have red squares surrounding the population number. SO-171 is a putative 
hybrid population between G. howellii and G. squarrosa. Map created by Jamul Hahn. (b) All 
populations used in this study, including Blue Mountain (Grindelia squarrosa) and Clear Creek, 
ID (G. howellii). Map from Google Earth. 
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a. | b. 

Figure 2 — Principal Coordinate Analysis ordination of Nei’s genetic distance for (a) All 
populations (first two axes explains 65.95% of the variation) and (b) Montana G. howellii 
populations (first two axes explains 84.22% of the variation). 
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Figure 3 —- STRUCTURE results for the three datasets. TOP: All 11 sampled populations 
including Blue Mountain (Grindelia squarrosa) and SO-171 (putative hybrid of G. howellii and 
G. squarrosa) (K=4). MIDDLE: Grindelia howellii populations, which excludes the putative 
hybrid SO-171 and Blue Mountain G. squarrosa populations (K=2). BOTTOM: Montana G. 
howellii populations (which excludes SO-171, Clear Creek, and G. squarrosa populations; K=3). 
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